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Abstract 

We consider the stationary one dimensional Schrodinger-Poisson system on a bounded 
interval with a background potential describing a quantum well. Using a partition function 
which forces the particles to remain in the quantum well, the limit h — y in the nonlinear 
system leads to a uniquely solved nonlinear problem with concentrated particle density. It 
allows to conclude about the convergence of the solution. 
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1 Introduction 

The quantum state of a gas of charged particles is described, in the mean field approximation, 
by a nonlinear one-particle Schrodingcr equation where the electrostatic repulsion is modeled by 
a nonlinear potential term depending on the charge density through a Poisson equation. This 
class of models is usually referred to as Schrodinger-Poisson systems. In this work we consider a 
stationary Schrodinger-Poisson system in a bounded region of R, for which a background potential 
models a quantum well, while the nonlinear potential extends on a wider scale. After introducing a 
rescaling for which the small parameter h > represents an inverse length scale, the support of the 
potential well squeezes asymptotically to a single point in the limit h — > 0. An equilibrium state 
of a gas of charged particles confined in the quantum well will be considered, while the nonlinear 
electrostatic potential created by such a concentrated charge tends to a potential which is picewise 
linear and almost constant in the well. 

Such a Schrodinger-Poisson problem has recently been considered in [3] , [4] and [14] in a more 
complex setting involving far from equilibrium steady states. This one-dimensional analysis leads 
to a reduced model which happens to be very efficient in the numerical simulation of the electronic 
transport through semiconductor heterostructures, like resonant tunneling diodes, see [T] and [2]. 
Nevertheless, in the present work the Schrodinger-Poisson system has a unique solution and the 
analysis shows that the unicity is asymptotically preserved. It is not the case in [3], [3] and [T4] 
where hysteresis phenomena are predicted (e.g. in [11] and [TB] ). 

For the sake of simplicity we shall use a low energy-filter in the definition of the partition 
function / (see equation (|1.5[) below), that is the quantum states with an energy larger than the 
threshold es are not occupied. With such an assumption only the quantum states confined in 
the well have an effect on the nonlincarity. This is an important point to get asymptotically the 
macroscopic quantities. 

The semi-classical analysis of such a model was performed in [B] in dimension d = 2 and 3 and 
it appears that in the limit, the potential vanishes almost everywhere and produces a non null 
spectral perturbation. This apparent contradiction is solved through a rescaling with parameter 
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h and the potential behaviour depends on the Green function of the Laplace operator. For the 
ID problem, the nonlinear effect remains visible at the macroscopic scale in the limit h — > 0, this 
provides a non trivial approximation of the solution of the Schrodinger-Poisson system. 

Like in [3] , [I] , [S] and |14j the analysis will be a mixture of nonlinear apriori estimates combined 
with accurate semiclassical and spectral techniques (we refer to: [5], [5], [S] and [H]) adapted for 
potentials with limited regularity The outline of this analysis is the following. We end this section 
by introducing the model and by stating our results. In Section [2 we give some asymptotics 
for the spectrum of the linear Hamiltonian. In Section UJ we present preliminary results for the 
Schrodinger-Poisson system. The limit of the Schrodinger-Poisson system is done in Section @] 
The semi-classical analysis tools, necessary for the asymptotics, are given in Appendix lAl and [Bl 

1.1 The model 

Let f2 = (0, L) be an open bounded interval and U a non positive function in C^K) supported in 
the ball of radius one centered in the origin. 

For Xq £ fi, we define the potential with center xq and radius of order h > 



U h (x) = U (j-^j , zeO. 



Our analysis being concerned with the limit h — > 0, we can choose, without loss of generality, h 
small enough so that the support of U h is included in Q. In particular, defining with u> h the support 
of U h , we assume that uj h C n for all values of h below a suitable positive constant: h < ho. 

Next we assign the function / <S C°°(R), with a threshold at £5 < and fulfilling the conditions 

f(x)>0, Wx<s s , (1.1) 

f(x)=0, Vx>e s , (1.2) 
f'(x)<0, VxeR, (1.3) 
and address, for h £ (0, ho], the following problem: find V h solving the nonlinear Poisson equation 

-*v»-.n mo (14) 

v "\m = ° 
where the source term is 

n[V h }=J2M)\^\ 2 , (1-5) 

i>l 

and {si/^j are the eigenvalues of the nonlinear Hamiltonian 

H h = -h 2 ^ + U h + V h (1.6) 
dx z 

numerated from inf a(H h ) counting multiplicities, while {'I'i 1 }^ ar e the corresponding eigenvec- 
tors 

{ ,, (1-7) 

W 1 = 

l Ian • 

The equations fj X .4[) . (jl.5p and (jl.7p define the stationary Schrodinger-Poisson system associated 
with the potential well U h and the function /. In practical applications, where these equations 
are used for the description of the charge distribution in electronic devices, n[V h ] describes the 
density of the charge careers of the system, while / is a response function which depends on the 
characteristics of the device and has to be considered as a data item of the problem. 
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The stationary states form a set of real normalized functions 

lm^ = 0;||^|| i2(n) =l, (1.8) 

The analysis of our Schrodinger-Poisson system, will involve the operator 

H = —^ + U; D(H )=H 2 (R), (1.9) 

whose point spectrum, a p (Ho), contains a finite number of points embedded in [— ||J7|| Loo ,0). In 
particular, we make the following assumption 

e x := infer (H ) <e s . (1.10) 

The hypothesis (|1.10[) - which prevents the solution to (|1.4[) - (|1.7|) to be trivial - will be extensively 
used in this work. 



1.2 Results 



In our one dimensional case (see [12]) and in the general case of the dimension d < 3 (see [13]), it 
has been proved that, for h > 0, the problem (|1.4[) - (|1.7[) admits an unique solution, V h in our 
notation, in Hq(Q). However, the uniqueness of V h is by far not obvious if the limit h —> is 
considered. The aim of our analysis is to understand the asymptotic behaviour of the system (|1.4D 

- (|1.7|) as h — > 0. This in order to provide a simplified modelling for the nonlinearities produced 
by charged particles confined in quantum layers. Such a program has been carried out in [6] in 
dimension d = 2 and 3. 

In our one dimensional framework, the functional spaces for the solution verify better Sobolev 
injections than in the 2D and 31? case, which gives, up to extraction, strong convergence results 
for the potential (see [3] and [4]). We have then enough regularity to obtain the convergence of 
the negative spectrum. This allows to determine the limit of the Shrodinger-Poisson system p. 41) 

- (|1.7[) . up to extraction, and the convergence of all the sequence follows from the uniqueness, and 
the explicit computation, of the solution of the limit problem. 

Our main results below, whose proofs are given in Scction[4j gather the asymptotic informations 
for the ID case. 

Theorem 1.1. The nonlinearity {y h )h£{o,h \ *s bounded in W 1,oo (0,L) and tends, strongly in 
C°' a (0, L), Va G (0, 1), to the potential Vq given by 



V {x) 



. i>l 



L 

(1.11) 



]T/(e !: + #) ] ^-(L-x), x <x<L 



, i>l 



L 



where {e,}i<i<Ar is the discrete spectrum of Hq, ei>N+i = and 6 £ (0,es — ei) is the unigue non 

negative solution of the nonlinear eguation 9 = Xq(1 — -j-) /(e.; + 9). 

i>i 

The density («[V r ' l ])^ g (o i h ] tends to the mesure 

M = E^ e * + ^°- ( L12 ) 

i>l 

for the weak* topology on the space A4b{0,L) of bounded mesures on (0, L). 



Hence for the ID problem, the nonlinear effect produced at the quantum scale remains visible 
at the macroscopic scale in the limit h — ► 0. 
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2 Some asymptotics for the linear operator 



Consider the operator Hq denned in fj 1 .9|) . The multiplication operator U is a relatively compact 

,2 

perturbation of — ■js- Then, the Kato-Reillich and Weyl Theorems imply that Hq is self-adjoint 
and a ess (H ) = [0, +oo). 

The reader may refer to Proposition 7.4 in [TB] to see that a p (Ho) ^ is always true in the 
one dimensional case when U < and not identically zero, and the discret spectrum of Hq is a 
countable set of negative eigenvalues with multiplicity one and zero as possible accumulation point. 
In addition, the Proposition 7.5 in [T8] gives the following estimation on the number N of negative 
eigenvalues of H : 

N < 1+ [ \x\\U(x)\dx, 

J Ft. 

bound which is finite in our case. It follows that 

a(H ) = {e u ...,e N }U[0,+<x>) (2- 1 ) 
where the ej are the points of ad(Ho). We will work with the convention 

e,i =0, for i > N + 1 (2.2) 
The present section is devoted to the description, when h — > 0, of the spectrum of the operator 

H' a l = -h 2 ^ + U h , D(H%) = H 2 r\Hl(Q) (2.3) 

using a comparison with the spectrum of Hq. The operator Hq is the linear part of the Hamiltonian 
H h involved in the Schrddinger-Poisson system (|1.4[) - (|1.7[) and, therefore, it will play an important 
role in the analysis: the spectral asymptotics of the linear operator will give information about the 
nonlinear one. 

We will denote by (— /i 2 ^)n the realisation of —h 2 -^ on fl with domain Hq n H 2 (Q). The 
spectrum of the operator (— h 2 -^)^ is the sequence of the = h l L 2 and the corresponding 

normalized eigenvectors are the (fi(x) = y^"sin ^ for i > 1. 
By setting 

A:=\\U\\ LOO > 0. 

we have: 

<r{{-h 2 ^-K)n) = {\ h i -Ki>l} 
Let e be a constant such that e > —A, we have: 

# (°((-h 2 ^ - A)n) n (-oo, s]\ = #{* > 1; - A < e} < VF+A^ 

Coming back to the operator Hq, we note that: 

and therefore Ve > — A, the integer 

N h l :=#(a( J ff fc )n(-oo, £ ]) (2.4) 

verifies 

^V' 1 < # (*((-h 2 -^ - A)n) n (-co, e]\ - 0^). (2.5) 

The previous asymptotic order is an important point in the proof of the following Lemma, which 
is the main result of this section. 
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Lemma 2.1. Consider e E (ej\r,0) and N h given by If we define e\ to be the eigenvalues of 

Hq, then, for ho small enough, we have: 

• N h = N, V/i € (0, ho] 



• Fori = 1, ...,N: 



lim e) 



Remark 2.2. A consequence of this Lemma is: 

liminf e< > 

for i > N + 1. In other words, the eigenvalues of Hq , with number i > JV + 1, are asymptotically 
embedded in the continuous spectrum of Hq. 

The proof of this lemma is based on the exponential decay of eigenfunctions given in Appendix 
iBl The estimates are written with weight functions involving Agmon distances which proprieties 
are recalled in Appendix O 

Proof. The states corresponding to energies below e are exponentialy decaying outside the support 
uj h of U h , therefore the behaviour of the eigenvalues of Hq is well discribed by the comparison 
with the operator considered on the whole space, Hq, defined as follows: 

Hq = -h 2 ^ + U h ; D(H' l ) = H 2 (R) (2.6) 
The operator Hq is unitary equivalent to the operator Hq through the unitary map on L 2 (R): 

4>{x) = h^(j)(hx + xq) 

and (|2.1|) implies that: 

a(H^) = {ei,...,ejv}U [0, +oo) 
We will denote by 4>f and $ the normalized eigenvectors of Hq and Hq respectively. 
• Introduce the family: u\ £ D(Hq) defined by 

u?= X #, fori = l,...,JV (2.7) 

where \ G Cq°(M) is such that: x = 1 on B(xq, -j) and x = on M\B(xo, R) for a raduis R > 
verifying B(xq,R) C ft- 

We have the following classical inequality for self-adjoint operators, which can be found in [8]: 

\\ut\\ LH n)d(ei,o-(HS)) < \\(H% - e^H^n) . 

We refer to [5] to see that estimation (|B.1[) can be extended to the whole space where the weight 
function is given by the Agmon distance on M related to the potential (U h — e). This, combined 
with inequality (|A.2[) for f2 = Ft, gives 



dx 



I" 

e C0 ~ 



i 

dx <C (2. 



and it follows that ||M^||L 2 (n) is bounded from below. Indeed, we have: 

h\\2 \lh\2j„ , k,|2|7/i|2 



K\\Un) = / m z dx+ / \xV\PiVdx 

JB(x ,R/2) JM\B(x ,R/2) 



h,2 dx + / (|x| 2 - Wl\ 2 dx = 1+ (|x| 2 - l)\$\ 2 dx 

I R\B{x ,R/2) J R\B(x ,R/2) 
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And, the decay estimate (|2.8[) implies that 



/ (|X| 2 ~ l)\^\ 2 dx < (Hxllloc + 1) / e -2c |s-*o|/fc| e co|*-*o|/^|2<fc 

I R\B(x ,R/2) J m\B(x ,R/2) 

< C [ \e Colx - Xol/h $\ 2 dx e - CoR/h < C'e- CoR/h 

which leads to 



IK h IU»(n) > J. V/iG(0,/i ] 
for ho small enough. Then, a direct computation gives: 

(fl* e,)«? = -h 2 X "% - 2h\'($y. 

The r.h.s. in the previous equality being supported in a region where ip? is exponentially decaying, 
the estimate (|2.8[) gives: 

\\(H£ - e.OulW^ <Ce~Z (2.9) 

where C, 7 > doesn't depend on i and h. 
We deduce that for i = 1, N 

d^a^)) — >0 
ft— >o 

It implies N h > N for all /i e (0, fr ] and h Q small. Then, it's enough to show that N h < N to 
conclude. 

• For i = 1, JV , if we set = x^ 1 £ D(Hq) using (|B.4[) . we can reproduce the calculation of 
the first point to obtain 

\\(H^-e^\\ LHM) <Ce-i (2.10) 
Now, let / = [— ||J7||i,=°,e] and a > be small enough to have 

<r(Hft) n ((/ + B(0, 2a))\I) = 0. 

We consider the vector space E spaned by v^, ...,trL h and the spectral subspace F corresponding 
to <j(Hq) n /. Estimation (|B.4[) implies that the matrix M = ((u'\ v , J l )i J 2^) 1<i _j <N h verifies 

M = I + 0{e~i) (2.11) 

when h — > 0. As it will be clarified in Remark 12.31 equation (|2.11[) implies that the v£ are linearly 
independant. If we consider in addition equation (|2.10p . the conditions for the application of 
Proposition 2.5 in [3] are gathered and the distance d(E,F) (definition given in [§]) is estimated 
by 



d(E,F)< 



a 



where \ m i n is the smallest eigenvalue of M . We deduce from (|2.1ip that A m j„ = 1 + o(l), and 
using (|2.5p . we get: 

< C(N h )i e -% < Ch~h-i < i 

for all /i € (0, ho] and /io small enough. This last condition allows us to state that the projection 
Hp\ E : E ^ F is injective (e.g. in Lemma 1.3 in 0), from which the condition N h < N follows. 
■ 

Remark 2.3. The matrix M := ((v^,Vj)) verifies M = J + (9(e h) for a constant 7 > where 
the asymptotics O is considered with respect to the norm 

||A||a = max \a i3 \ 
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This implies M € GL N h(C) and the linear independence of the v!?. Indeed, remarking that (see in 
®) 

VAeC NhxN \ ||A|| 2 <iV ,i ||A||A 

where 

xjtO N| 2 *-f 

we obtain ||M-/|| 2 < Cft -1 e~* <l,Vfee (0, ho] for ho small. The inequality ||AB|| 2 < ||A|| 2 ||B||2 
implies that the sequence S := Yl n >o(~^) n (M — I) n converges and verifies: 

SM = MS = 1 

and the matrix M is invertiblc. Now, suppose there exists A = (Aj)i such that \iV^ = 0, then 
X T M = and A = 0. 

The asymptotics \ m in = 1 + o(l), where A TO j„ is the smallest eigenvalue of M, is given by the 
perturbation estimate for eigenvalues 

|A„(M)-1| < ||M-/|| 2 , a(M) = {A„(M), 1 < n < N h } 

which can be found in 0. As ||M - /|| 2 < JV fc ||M - 7|| A < Ce~*, Vft G (0, h ], we get: 

Xmin = H-0(e~*) 



3 Preliminary results for the Schrodinger-Poisson system 

We give some apriori estimates for our Schrodinger-Poisson system l|1.4|) - (|1.7|l . The results 
presented below are also true in dimension d < 3 as it appears in [6]. 
From the lower bound 

d 2 

—r^V h > in fi 

with homogeneous Dirichlet boundary conditions, the maximum principle implies 

V h > in fl (3.1) 

Thus, V h defines a positive perturbation of the Hamiltonian Hq given by (|2.3[) . The spectra of 
is bounded from below by the norm H^H^^r) and we can state 

mf{4} i > 1 >-||f/|U • (3.2) 

Due to the definition of the source term (|1.5[) . y' 1 is generated by those energy levels e\ placed 
below the cut off ss of the characteristic function /. In order to study the semiclassical behaviour 
of our system, we are interested into the spectral properties of the Hamiltonian 

d 2 

H h = -h 2 — + U h + V h , D(H h ) = H 2 n H^{n) (3.3) 

ax 2, 

in the spectral interval [— ||C/|| ioo , £5), as h — > 0. In particular, a uniform bound for the number 
of eigenvalues e 1 } £ [— ||C/|| LOO , es) as h — > is required. As noticed above, the operator H h is 
obtained as a positive perturbation of Hq through the Poisson potential V h . Then, the minimax 
principle implies that 

e 1 } < £?, V^ > 1 (3.4) 
and the Lemma [2~T1 leads to the following result. 
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Lemma 3.1. Let N be the integer given by (|2.1[) . then Vh € (0, ho] 

#(a(HS)n[-\\U\\ Lao ,e s ))<N (3.5) 

#(<r(H h )n[-\\U\\ L „,e 8 )) <N (3.6) 
where o~(H) denotes the spectrum of H. 

We recall the variational formulation, given in |13| in dimension d < 3, of the Schrodingcr- 
Poisson problem (|1.4|) - (jl.7D . Rephrasing the results of this work for our system, we can state 
that the solution to the equation (jl.4|) - (jl.7|) is equivalent to the minimization problem 



inf J(V); J(V) = \( \V(x)\ 2 dx + Tr [F (H h (V))] , (3.7) 
where F is the positive function 

r + co 

F(x)= / f(s)ds, 

J x 

while the Hamiltonian H h (V) is given by 

H h (V) = ~h 2 ^ + U h + V, D{H h {V)) = H 2 n H^(fl) . 
dx z 

Moreover, the function J(V) is Frechet-C°° w.r.t. V, strictly convex and coercive on H^fl) and 
p.7p admits a unique solution in this space. The following Proposition is a direct consequence of 
this result. 

Proposition 3.2. The solution to the Schrodinger-Poisson problem - {1. 7\ ) is bounded in 

Hq(Q) uniformly with respect to h. 

Proof. From the variational formulation recalled above, the solution V h is the minimum of the 
convex map J{V), therefore we have 

f \{V h Y{x)\ 2 dx + Tr [F (H h {V h ))] < J(0) = Tr [F (H h (0))] 

where H h (V h ) simply coincides with the Hamiltonian H h , while H h (0) can be identified with Hq 
defined in (|2.3[) . The relation 

N 



Tr [F(H h (V h ))] =^F(4)>0, 



i=l 



with N given in Lemma 13. 11 implies 



V h \L m < 2T r[F(H>>)] . (3.8) 



From Lemma 1 2. 1[ the explicit expression of the r.h.s. here is 

N 

Tr [fK)]=^F(ef). 

i=l 

The result easily follows by combining (|3.8[) with the inequality 

N 

Y J F(e'l)<N sup F< 

i=l xe[-\\U\\ L ^,e s ) 



Next we use the assumption (|1.10[) and Lemma 12.11 to get uniform upper bound for the first 
spectral point of H h as h — s- 0. 
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Lemma 3.3. For ho small enough, the condition 

e\ < e s (3.9) 

holds for all h £ (0, ho] ■ 

Proof. We use a reductio ad absurdum argument. Let h £ (0, ho) be such that e\ > es- It follows 
from (|1.2[) and from the definition (|1.5[) that the corresponding charge density, n V h , and, then, 

the Poisson potential V h are null in il. In these conditions the Hamiltonians H h and Hq coincide 
and we have 

4 = A > es ■ (3.10) 

On the other hand, as it follows from Lemma [2.11 we have e\ — > inf a (Hq) when h — > 0. Then 
from the assumption (|1.10[) . the condition 

e\ < e s (3.11) 
definitely holds for fc->fl, which is in contradiction with (|3.10[) . ■ 



4 The semi-classical limit 

In this section, the asymptotic behaviour of the Schrodingcr-Poisson system (|1.4|) - (|1.7p is de- 
scribed. Contrary to the other sections, the results presented here can not be extended to higher 
dimensions. Therefore, we will use the notation (0,L) introduced in Section Q] to design the open 
domain O. 

Due to (|1.2j) . only a finite number of eigenvectors contributes to the density and therefore n[V^ /l ] = 
/(£^)|\&^| 2 belongs to £ 1 (0, L) as a finite sum of L 1 functions. Then, taking into account the 

i>l 

embedding i? 1 (0,L) C C(0, L) in dimension 1, the accurate function space for the potential is: 

W 2 (0, L) := {V £ C(0, L); V" £ M b (0, L), V(0) = V(L) = 0} 
The space BVq(0, L) with the norm: 

\\V\\bv :=||V|| o(0)i ) + |ini« 

is a Banach space. Here ||^|| m — |>Lt| (0, -L) is the strong norm of the bounded mesure fi which is 
equal to the L 1 norm when [i £ L 1 (0, L). 

We have the following continuous injections which will be useful in this work (e.g. in [3], [3]): 
Va G (0,1) 

W 2 (0,L) W l '°°(0,L) t W 2 (0,i) CC C°' a (0,L) (4.1) 

where the second injection is compact. The continuous injection BVj^O, L) c — ^ W 1,o °(0, L) doesn't 
present real difficulties, the distribution function of a bounded mesure being regular. The second 
injection is a consequence of the compact embedding W 1,oo (0,i) CC C°' a (0,L). 
We can already give some apriori estimates for the potential and the density. 

Proposition 4.1. The density {n\V h ])h^(Q ^ is bounded in Aib(0,L). The potential (^' i )/ l e(o./i ] 
is bounded in W 1 ' 00 ^, L) and relatively compact in C 0,a (0, L) for any a £ (0, 1). 

Proof. We start looking for an estimate on the density. Using the normalization condition (|1.8[) , 
we have: 

\\n[V h ]\\ m = [ L ' n[V h }dx = J2f(£iWi\\Uo,L) 

J o i>l 
< N h sup f(x) 
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where N h denotes the number of eigenvalues of H h below £5. According to Lemma T3. 11 we have 
N h = 0(1) when h tends to and the density is bounded in Mb{0, L). 

For what concerns the potential V h , the Hq bound provided by Proposition ^. 2l and the continuous 
embedding Hq(0,L) ^ C(0,L) implies 

\\V h \\c(p,L) < C, Vhe(0,h o ] 

Then, it follows from the Poisson equation — (V )" = n[V h ] and the bound on the density that 

\\(V h )"\\ m = \\n[V h ]\\ m <C, Vhe(0,h o ] 

We obtain that (V )he(o,h ] * s bounded in BV 2 (0, L) and (|4. 1[) allows to conclude. ■ 

The density {n[V h ]) h€ ^o,h Q ] being bounded in Mb(0,L), it is relatively compact in the weak* 
topology. Therefore, we deduce from Proposition 14.11 that for any < a < 1 , it is possible to 
extract from any infinite set S C (0, ho], which has as an accumulation point, a countable subset 
D such that £ D and: 

]im(n[V h ]-fx, 9 )=0, VtpeC([0,L]) (4.2) 

h— >0 
h£D 

Um\\V h -V o \\o, a = 0, (4.3) 

h->-0 
heD 

for some n € Mb{0,L) and Vq G C°' a (0,L). Here, ||.||o,a is the usual norm on the Holder space 
C°' a (0, L). Then we have the following result on the convergence of the spectrum of H h : 

Lemma 4.2. Consider < a < 1, D C (0, ho] s.t e D, and V € C°- a (0,L) verifying ffty . If 
9 = Vo(xo) and N\ is the greatest integer such that ejVj + < 0, then: 

• 6 > and N-l > 1 

• For i = 1, iVi 

lim e' 1 = Si + 9 

h^>0 
heD 



• For i > Ni + 1 



liminf £ > 

heD 



The proof of this lemma is similar to the proof of Lemma 12.11 where the nonlinearity have to 
be considered as it is done in |15) and [5] in the two dimensional case. 

Proof. The asymptotics of the ej 1 is given by the comparison with the Hamiltonian: 

He =-h 2 ^ + U h + 9, D(H' g l ) = H 2 (M) 

The operator H@ is unitary equivalent to Hg defined by: 

d 2 

H e = -— + U + 9, D(H e )=H 2 (R) 

and therefore: 

a(H> 1 ) =K,...,4}U[0,+co[ 

where e\ := ej + 9. We note that the normalized eigenvectors of Hg are the cj>^ introduced in the 
proof of Lemma 12.11 

The W 1 ' 00 bound obtained in Proposition 14 . II implies that there exists a constant C > such that 
V/i 6 (0,ho]: 

|^ h (i)-V fc (|/)|<C|a;-|/|, Vx,yen (4.4) 
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and the condition (|4.3[) gives: 

lira \V h (x ) - 6\ = (4.5) 

h— >0 
h£D 

The previous equations are important in our nonlinear framework. Indeed, they asymptotically 
enable to consider V h as a constant on a domain outside which the modes of interest are exponen- 
tially decaying. 

Note that 9 > is a direct consequence of (|3.1j) and (|4.5[) . 
• Set e e (-||C/|| L oo,0) and define: 

iV /l = #(a(i/ h )n(-oo, £ ]) (4.6) 

We introduce the family: v\ = € D(H' g l ), for i = 1, TV' 1 , where \ G Co° W is s -t- X = 1 on 
B(aio, "f ) an d X = 011 K\S(xo, i?) for a radius R > verifying B(xq,R) C fi. 
As seen in the proof of Lemma 12.11 the exponential decay (|B.4|) for implies that it is enough 
to estimate (Hg — e^)v^ to get 

d(e?,a(H$))l i<m — > (4.7) 

n— >-0, ft£D 

Now we have 

(H> 1 - £ >>^ - -h\"^ - 2h 2 x '(^)' - x (V h - 

The function x being regular with derivatives supported where the eigenfunction is exponentially 
decaying, we have: 

\\h 2 X '^1 + 2h 2 X '(^y\\ LHM) <Ce-* 

for C, 7 > 0. Next, we define R h = ft-hi^ and we suppose that ho is small enough so that 
B(xq, R h ) C £1. Proposition 14 . 1 1 implies that the potential V h is bounded in L°°(SY), therefore, the 
application of the estimate (|B.4I) leads to: 



\x(v h -mn"<\\ x \\i 



If \(V h -9)^\ 2 dx+ f \(V h -9)^\ 2 dx 

\JB(x ,R h ) Ju\B(x ,R h ) j 



<c[\\V h -9\\l^ B(xatRh)) + jje*^^\ 2 dxe- 2 ^ 

<c\\v h -e\\l x{B{xo>Rh)) + o{e-*^) 

On the other hand, we deduce from equations (|4.4p and (|4.5[) that Vx <E B(xo,R h ) 

\V h (x)-9\ < \V h (x)-V h (x )\ + \V h (x )-9\ <Ch\n\+o{l) 

h 

when h — > 0, h £ D. Here, the asymptotics o(l) doesn't depend on x and therefore 

lim\\x(V h -9)^\\ L2m W =0. 

heD 

We remark then that a function a{h) independant of i can be found such that a(h) — > when 
h — > and 

IK^-e^lUafRjW* <a(/i), VfceD. (4.8) 



and <|477j) follows. 

According to Lemma \'3.'Sl we have e\ < £s, therefore the application of (|4.7|) with e = £5 implies 
< and iVj. > 1. This provides the first point of the Lemma. 
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• Consider the family it' 1 £ D{H h ) given by (|2.7[) . Then, for i = 1, ...,Ni, we have: 

(if* - ef)«? - -/»V# - + X (U' 1 - 0)# 

Using the exponential decay estimate (|2.8[) for (j>^ , the same computations as in the previous point 
can be performed to obtain 

d(ela(H h )) -> (4.9) 

ft— >0, h£D 

In what follows, we set e £ (e^ , 0) and N h the corresponding integer given by relation 
From Lemma 12.11 we have N h = 0(1). Then, taking into account the results (|4.7[) and 
above, we can already deduce that the following statement holds: for all 5 > small, 3ho s.t. 
\/h £ D n (0, /io] we have, even if the numerotation of the eigenvalues has to be changed, 

4,i,-,4,N k £l4-S,4 + S] for k = l,...,N 1 (4.10) 

where Nk > 1 is an integer depending on h such that = N h . This implies N h > N\ and 

if we obtain N h < Nx, Vh £ D PI (0, /io], the proof will be completed. 

• Consider I = [— \\U\\loo, e] and a > small enough so that 

<r(H$) n ((/ + B(0, 2a))\I) = 0. 

Let E be the vector space spanned by v^, and F the spectral subspace associated to a(Hg)C) 
I. From the estimate (|B.4|) . the matrix M = ((^f , t'j)L 2 (.ZR))i<i,j<jV' 1 verifies 

M = I + <D(e-%) 

when h — >• and the are linearly independant. Then, equation (|4.8[) and Proposition 2.5 in [5] 
give 

where X m in is the smallest eigenvalue of M . As we already noticed, N h = 0(1) and therefore: 

d(E,F) < Ca(h) < - 

for all h £ D fl (0, /iq] with /iq small enough. Then, Lemma 1.3 in [3] gives N h < N-y. m 



In what follows, we give, for a set D verifying (|4.2p and (|4.3[) . the limit problem when h — > 0. 
Then the unicity of the limit allows to determine the asymptotics of the solution of (|1.4[) - (|1.7I) 
when h — > 0. 

Proposition 4.3. Consider < a < 1, L> C (0, /i ] s.t e I? 7 fi £ M b (0,L) and V £ C°' Q (0,L) 
verifying \4- r ity and |^73p. Then, we have 

^ = Y,f^ + V ^ x o)) S -o (4-11) 
i>i 

and the potential Vq is solution of the problem 

d 2 

~d^ V ° = ^ ( °' L) (4.12) 

Vb(0) = V (L) = 



Proof. Set = Vb(a;o) and take e £ (ejvi + 0) where A^i is the integer defined in Lemma 143 
Then, from Lemma |4.2[ we have for ho small enough: 

#{a{H h )C\ {-oo,s}) =N U Vh£Dn(0,h o ) 
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and the density can be written as follows: 



i=l 



/? I 2 



where e' 1 < e. Using the normalization (|1.8p and the exponential decay estimates, we get that 
| , 5' i | 2 (£E is of mass 1 and concentrates around xq. It implies that for i = 1, iVi 



l*' l | 2 dx 



(4.13) 



for the weak* topology on the space of bounded mesures on (0, L). Indeed, for i = 1, ...,N\ we can 
apply the estimate (|B.4[) to write 



/V'f Vdx 


< f \<p\ e" 2co|a 


;—x \/h 


e c \x-x \/h-qjh 


2 

dx 


Jo 


J supp <p 









< \\f\\ 



D c \x-x \/h-qjh 



sup e -^o\x-x Q \/h < c\\ip\\ Laa sup e 

L 2 (0,L) x^supp^ zEsuppip 



2c \x — xo\/h 



for any ip € L°°(0,L) and V/i £ flfl (0, ho]. It implies that, when supptp is a compact set in 
(0, L)\ {xq}, there exists c v > such that 



By taking the limit as /i — > 0, we get 



< C \\(p\\ LO o e- c ^ h , yheDf](0,h o ] 



hm 

h->0 
heD 







for all continuous function <p € C'QO, L]) with supp ip C (0, i) \ {a;o}- 
We recall that Lemma |4~21 also gives for i = 1, N%: 

lim £*! — + 

Therefore, from the unicity of the limit (|4.2p . we get the result (|4.1ip using the convergence f|4. 13[) . 
the continuity of / and the convention (|2.2I) . 

On the other hand, the convergence (|4.2[) , (|4.3p is valid in X>'(0, L): 

lim(n[^' 1 ] -/i,<^) =0, lim(^' i -F ,(p) -0, Wtp e C °°(0, L) 
ft,ez) dec 

From the Poisson equation — 4^V h — n\V h ], V7i G (0,/io] an d the continuity of the derivative on 
V'(0,L), we get 

-^o = M, 2^(0, L) 

As a consequence of (|4.3p . the potential V h tends to Vq strongly in C°' Q (0,L) and the boundary 
conditions appearing in the problem (|4.12p follow from 

Voix) = lim V h (x) 

heD 

for x = and x = L. m 

We conclude this section giving the proof of Theorem 11,11 
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Proof of Theorem 11.11 First, we recall that the IV 1,oo (0,L) bound for the potential was 
obtained in Proposition ^. II 

Next, let a be a constant in (0,1) and consider D C (0, h n ] s.t. <G D, [i <G Mb{0, L) and 
Vo £ C°' a (0,L) verifying (|4.2p and (|4.3p . As a consequence of Proposition 14. 31 the potential Vb 
verifies the problem (|4.1ip(|4.12p . and it can be computed explicitly as a function of 9 := Vo(xo) by 
solving the following transmission problem: 

f -(Vb)" = 0, (0,x o )U(x o ,L) 

\ Vo(0) = V (L) = 

with the jump conditions: 

r K,(xo) = y (4) 

) K(x+) - Vt(xo) = -Q2 f(ei + 0)) 

\ i>l 

The formula (jl.lll) for Vb follows. Then, equation (|l.ll|l considered at x = xq implies that the 
value 9 of the potential Vq at Xo, is solution of: 

e- Xo (l-^-)J2f(ei + 0) = (4.14) 

i>i 

According to the first point in Lemma 14.21 we have Vq{xq) > and we are interested in non 
negative solutions of (|4.14j) . If we call G(9) the l.h.s. of this equation, we remark that G is a 
continuous, strictly increasing function on [0, +oo) and such that lira G(9) = +oo. Therefore, it 

8^+00 

defines a bijection from [0, +oo) to [G(0), +oo). We have G(0) = -x (l - 7 1 ) J2i>i f( e i) < an< i 
G(es — ei) = £s — e\ > by the assumption (|1.10p . We deduce that there exists a unique value 
9 > solving (|4.14p . moreover, it verifies 9 £ (0,es — ei). 

It follows that the function Vq (resp. fi) verifying f)4.3[) (resp. (|4.2p ) is given in an unique way by 
(II. lip (resp. (|1.12|l ) where 9 is the positive solution of (|4.14l) . This gives the convergence results 
anounced in our theorem. 

Indeed, suppose that the convergence doesn't occure. Then, there exists a function (po £ C([0, L]), 
a constant e > and a set S C (0, ho] s.t. € S verifying: 

H^-Volka + K^' 1 ]-^)! >e, Vfces (4.15) 

By applying Proposition 14.11 we can extract a set D C S s.t. G D and (|4.2I) . f)4.3[) are verified 
for some functions fl € Mb(0,L) and Vq € C°' a (0, L). Then, according to the previous uniqueness 
result fl = /i and Vq = Vo and we get a contradiction comparing (|4.2p . (|4.3p with (|4.15p . ■ 

Remark 4.4. In the proof of Theorem 11.11 above, we have shown that the solution 9 of (|4. 14[) 
is in the interval (0,£s ~ ei). The condition 9 > implies that the sum X)i>i /( e « + ^) i s finite. 
Although the bound 9 < es — e.\ has no impact on the proof, it gives that ^2 i>± f{&i + 9) > 0, 
and therefore, that the limit potential Vb and density /i are not trivial. We deduce that Theorem 
11.11 provides a non trivial approximation of the solution V h of the problem (|1 .4[) - (| 1 . 7[) in the 
semi-classical limit h — > 0. As noted in the introduction, it is not the case anymore in dimension 
d = 2 and 3 as it appears in [6]. 
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completion of this study. He wishes to thank the support by the project "QUATRAIN" (No. 
BLAN07-2 219888, founded by the French National Agency for the Research). 
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A Agmon distance 



We will mainly follow [10] and present an additional result in the framework of our problem. Let 
/ be a real-valued function, continuous on a bounded connected set fl C 1R. 
For given x, y £ fi, the Agmon distance related to / is: 

d(x,y)= [ (fh(s)))lW(s)\ds 
Jo 

where 7 is the segment linking x and y. 

We remark that this distance is degenerated: it may happen that d(x, y) = with x ^ y. However 
the Agmon distance verifies the following properties: Vs, y, z £ ft 

d(x, y) = d(y, x), d(x, z) < d(x, y) + d(y, z) 

For y £ Q fixed, the fonction x <-> d(x, y) is Lipschitzian continuous. Then, it is differentiable 
almost everywhere by Rademacher's Theorem and everywhere x 1— > d[x, y) is differentiable, we 
have: 

\V x d(x,y)\<(f(x))l 

For a set w C ft, we define: 

d(x,ui) = inf d(x,y) 
Then, d(x,ui) has the same regularity as x <— > d(x,y) and 

\Vd(x,(j)\ < (f(x))l a.e.Q (A.l) 

In our case, we will use the Agmon distance associated to the potential U h — e for a constant e < 
which corresponds to taking / = U — e. 

In the exponential decay estimates, it will be useful to replace the Agmon distance with the 
Euclidian distance. To do this, the main point is to remark that for / = 1 the Agmon distance, 
d(x,y), corresponds to the Euclidian one, \x — y\. 

Lemma A.l. Let uj h be the support of the well U h , then 3cq,ci > such that: 

d(x,oj h ) > cq\x - xq\ — c\h, VieO (A.2) 

Proof. Let x G Q \ u> , y £ u h and consider the segment 7 : [0, 1] — s- ft linking x and y. Then, 
there exists a unique value so £ (0, 1] s.t. 7(so) £ duj h . Defining y = j(sq), we have: 

(tf h (7(*))-e)Il7'(*)l<k > / (U^^-eyiWWds 



= |e|a / \ 7 '(s)\ds = \s\*\y -x\ 
Jo 

Recalling that u h C B(xq, h), it follows: 

d(x,y) > |e|3(|a;-a;o| - h) 



B Exponential decay of eigenfunctions 

The exponential decay estimates, also called Agmon estimates, form a standard technical tool in 
evaluating the rate of decay of eigenfunctions far from the interaction support. In what follows, 
we apply this technique to the case of the Schrodingcr Poisson system with a squeezing quantum 
well; in particular we give some useful decay estimates for those stationary states related to the 
energies below some negative energy. 
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Lemma B.l (Agmon estimates). Let e h be a spectral point of the Hamiltonian 

H h = -h 2 A + U h + V h 

placed below some negative energy: e h < e where e € (— ||E/||l°°,0). The related normalized 
eigenvector admits the estimate 





+ 


e <p/h^h 




L 2 (n) 





(B.l) 



L 2 (0) 

where C is a suitable positive constant, is the weight junction 

<p{x) = {I- 5)d{x,u h ), xen 

6 is a positive parameter smaller than 1, d(x,y) is the Agmon distance introduced in AvvendixVfi[ 
while uj h is the support of U h . 

Proof. We use the relation (see for instance Theorem 1.1 in [§]) 

h 2 [ \\7(eiu)\ 2 dx+ [ {V - \\7<f>\ 2 )e 2 i u 2 dx = f e 2 i (-h 2 A + V)u-udx 
Jn Jn Jn 

Setting u = 4> = (1 - S) d(x, w h ) and V = U h + V h - e h , we get 

h 2 f |V(e*¥ h )| 2 dar + f (U h + V h - e h — |V0| 2 )e 2 ^ {V h ) 2 dx = 
Jn Jn 

Next we follow the same line as in Proposition 3.3.1 of [8] and introduce the set 

n+ = {xe{l\U h -e>5} , 

flj = {x e Q I U h -e < 6} , 

where 6 is a positive parameter such that 6 < min(l, \e\). As it appears in (|3.ip . we have V h > 
and therefore 



h 2 f \V(e^ h )\ 2 dx + f (U h + V h - e h - |V0| 2 )e 2 £ {^ h ) 2 dx = / 
Jn Jn+ Jn\ 



(|V0| 2 - {U h + V h - e h ))e 2 i (V h ) 2 dx 



< I (|V^| 2 - (U h -e))e 2 i {^ h ) 2 dx < f ((U h - e)+ - {U h - e))e 2 i {^ h ) 2 dx 

In — In — 



(U h - e)_e 2 ^ (^ h ) 2 dx < 



L oo / e 2 £(V h ) 2 dx 
nj 



where we used e h < e and the inequality 

\\74>\ 2 <(i-5) 2 (u h -s) + , p.p.n 



(B.2) 



(B.3) 



which follows from (jA.ip . 

On Clg, we have (U h — e)+ = U h — e. Using (|B.3j) again, we obtain that \tx £ Clf 
U h + V h - e h - |V</>| 2 > U h - e - |V0| 2 > U h - e - (1 - 5) 2 {U h - e] 
= S(2 - 8)(U h - e) >S 2 
Injecting the result above in (|B.2[) . we get: 

\\hV{ei^ h )\\ 2 LHa) +8 2 [ e 2 i {^ h ) 2 dx < \\U\\ L ~ f e 2 % (V h ) 2 dx 

I n+ In — 
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The domain £l g is inside the support of U h ; therefore, in this region, we have <j> = and taking 
into account the normalization condition H^'Hl^o) = 1, it follows: 



||W(e**' l )||i 2(f2) +^y e^(^ l Ydx< (||[%°= +5 A ) \ e^ ^> h ydx 
This gives the estimate (|B.1[) . ■ 



< 



Corollary B.2. Under the assumptions of Lemma (|B.1[) . the following estimate holds 



I x ~ x O I i. 



+ 


1 x — x O 1 L 


L 2 (f2) 





L 2 (n) 



(B.4) 



/or suitable positive constants Co and C. 



Proof. As a direct consequence of the estimate (|B.1[) and the inequality (|A.2|) . it follows that the 
L -norm of the function e c ° s $ is uniformly bounded w.r.t. ft. 
For what concerns the first term in (IB.4I) . we notice that 



|V*' l | < 



The term Ve^ at the r.h.s. is pointwisc bounded by 

1 



< 



£ )'+ eh 



as it comes from (|A.1[) . Then, using once more the relation (|B.1[) . we obtain 



e" Cl 


fte c ° k V*' 1 


< 








L 2 (f2) 





+ sup (U h -e)l 



L 2 (Q) 



<c. 
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